Disappearance of integer quantum Hall effect 
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The disappearance of integer quantum Hall effect (IQHE) at strong disorder and weak magnetic 
field is studied in a lattice model. A generic sequence by which the IQHE plateaus disappear is 
revealed: higher IQHE plateaus always vanish earlier than lower ones, and extended levels between 
those plateaus do not float up in energy but keep merging together after the destruction of plateaus. 
All of these features remain to be true in the weak-field limit as shown by the thermodynamic- 
localization-length calculation. Topological characterization in terms of Chern integers provides a 
simple physical explanation and suggests a qualitative difference between the lattice and continuum 
models. 
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It is an important issue how the integer quantum Hall 
effect vanishes in the weak magnetic field limit. Several 
years ago, Khmel'nitzkii Jlj] and Laughlin D] both ar- 
gued that extended states at the centers of Landau levels 
would not disappear discontinuously nor merge together, 
but should rather float up indefinitely towards higher en- 
ergy to pass the Fermi level. This argument is also a key 
scenario in the global phase diagram proposed by Kivel- 
son, Lee and Zhang ||] for the quantum Hall effect. As 
a consequence, a direct transition from a higher IQHE 
state to an insulator state is forbiden. || Recently such 
a floating picture for extended states has been challenged 
by a numerical study by Liu et al. [|J based on a tight 
binding model (TBM), in which extended levels as rep- 
resented by peaks of a finite-size localization length are 
found to disappear without changing their position at 
strong disorder. It indicates that the transition of IQHE 
states to an insulator occurs without involving a float- 
ing up of extended states. However, Yang and Bhatt || 
recently argued that this feature in a lattice model may 
not be contradictory to the conventional floating picture 
in the continuum limit. The argument goes as follows. 
They observed that in the TBM the critical strength W c 
of disorder by which all extended states vanish is roughly 
W c ~ 6t (t is the hopping integral) at different magnetic 
flux strengths. If such a finite W c indeed persists into 
the weak magnetic-field weaker disorder strength 

(W < W c ) will not be able to destroy extended states 
and the only way to get into an insulator phase in the 
weak-field limit would be a floating up of extended levels 
as suggested by Khmel'nitzkii and Laughlin. 

As these numerical calculations based on the TBM 
are all performed in relatively strong magnetic fields, 
one needs more conclusive evidences in order to resolve 
the above controversy and make the case relevant to the 
weak-field limit, which is of experimental interest. In 
this paper, we reexamine this problem in both strong 
and weak magnetic fields and present convincing results 
with regard to the physical mechanism of the disappear- 
ance of the IQHE and its behavior in the weak-field limit. 
Firstly, we show a systematic destruction of the IQHE at 



strong disorder, in which the IQHE plateaus are found 
to be destroyed in a one-by-one order from high to low 
energies. During such a process, extended levels sepa- 
rating different IQHE plateaus merge together and then 
disappear without floating up in energy, and the lowest 
extended level is the last to vanish at a critical disorder 
strength W c . Such a picture for destroying the IQHE 
is shown to persist into the weak field limit, where W c is 
found to decrease continuously with the decrease of mag- 
netic fields, and can be extrapolated to zero at zero field 
limit, contrary to a field-independent W c as suggested in 
Refs. |,§. A topological reason for such a destruction 
of the IQHE by disorder is given based on the Chern 
number which characterizes extended states: Q nonzero 
Chern integers with opposite sign are moving down from 
the band center to annihilate those in the lower-energy 
extended levels, and eventually a total annihilation of 
Chern numbers leads to a global insulating phase. In 
contrast to the finite-size localization length calculation 
in Ref. ||, though, we find that neighboring extended 
levels do merge together before their disappearance, as 
the result of moving and annihilation of Chern numbers. 
The TBM Hamiltonian is given as follows: 



H = 



<ij> 



'etc, 



H.c. 



E 



WiC Ci, 



where the hopping integral t is taken as the unit, and cf 
is a fermionic creation operator with < ij > referring to 
two nearest neighboring sites. A uniform magnetic flux 
per plaquette is given as <f> = J2a a v = 27r/M, where 
the summation runs over four links around a plaquette. 
And Wi is a random potential with strength < W/2 
(note that the disorder strength W here is twice bigger 
than that defined in Ref. || ) , and the white noise limit is 
considered with no correlation among different sites for 

Wi. 

An overall picture for the Hall conductance calcu- 
lated by the Kubo formula is shown in Fig. 1 with 
the flux strength <j> — 2tt/8 (only E < part is shown 
here). At weaker disorder (W — 1), three well-defined 
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IQHE plateaus at an = ve 2 jh {y — 1,2,3) are clearly 
shown, corresponding to four Landau levels centered at 
the jumps of the Hall conductance at E < 0. With W 
increasing from 1 to 6, one sees a systematic destruc- 
tion of these IQHE plateaus. For instance, at W = 2 
the v = 3 plateau is the first beginning to disappear, 
while at W = 3 and 4 the v = 2 IQHE plateau starts 
to vanish too. Meantime the lowest plateau still remains 
well defined at W = 4 and, in the inset of Fig. 1, such a 
v — 1 IQHE plateau region is partially shown at different 
sample sizes, where all the curves cross at a fixed-point 
corresponding to the Hall conductance \e 2 jh which is to 
be extrapolated to a sharp step of jump in the thermo- 
dynamic limit. By contrast, in those higher energy re- 
gions without the plateaus the Hall conductance is seen 
to monotonically decrease at larger lattice sizes. A break- 
down of the lowest IQHE plateau is eventually found at 
larger W's (W = 6 case is shown in Fig. 1). Hall con- 
ductances at weaker flux strengths with M = 11, 16 and 
24 all exhibit the similar generic features as in Fig. 1. 

A key thing to understand the above evolution of IQHE 
plateaus is the localization-delocalization transition. It 
is a well-known fact that in the IQHE an extended state 
can be characterized by a topological index, namely, a 
non-zero Chern integer. [Q And the boundary-condition- 
averaged Hall conductance is a summation of all those 
Chern numbers carried by states below Fermi surface: 
§ < cr H {E f ) >= e 2 /hJ2 em<Ef C (m) > where the Chern 
number is always an integer. Since its distribu- 

tion in the thermodynamic limit decides delocalization 
regions, one may define states with nonzero Chern in- 
tegers as extended states and calculate the correspond- 
ing density of extended states p ex t- || The results are 
presented in Fig. 2 with the same flux strength as in 
Fig. 1. At W = 1 (weak disorder), well-defined peaks 
of p ext are located at centers of Landau-level bands, sep- 
arated by localized regions represented by plateaus in 
Fig. 1. Widths of these peaks will approach to zero 
in the thermodynamic limit, [pi Total Chern number for 
each of three lower-energy peaks is found to be exactly 
+1, which is the reason leading to three quantized Hall 
plateaus at +1, +2, and +3 in unit of e 2 /h shown in 
Fig. 1, when the Fermi surface is located between these 
peaks. The last peak closest to the band center carries a 
total Chern number —3, which guarantees that the Hall 
conductance in Fig. 1 falls back to zero beyond the third 
plateau when the Fermi energy approaches E = 0. This 
is a peculiar feature of a lattice model since when the 
whole band is half-filled, the Hall conductance has to be 
zero. These Chern integers are conserved as long as they 
are well separated by localization regions, just like the 
quantization of the plateaus themselves. 

The distribution of Chern integers (thus the extended 
states) begins to qualitatively change at stronger dis- 
orders. At W = 2 and 3, two higher-energy peaks in 



Fig. 2 start to merge due to disorder scattering. Be- 
cause these two peaks carry total Chern numbers of +1 
and —3, respectively, a mixing and annihilation of these 
Chern numbers results in a disappearance of the third 
plateau in Fig. 1, and leads to a substantial reduction 
of the magnitude of the Hall conductance in that region. 
By contrast, the lowest-energy peak remains separated 
from the rest spectrum up to W = 5. In the insert of 
Fig. 2, such a peak at W — 4 as located at the energy 
E c i ~ —3.4 is shown to become narrower and sharper 
with the increase of lattice sizes, which indicates that a 
well-defined mobility edge still exists at E c i . Such a peak 
disappears finally as the negative Chern numbers reach 
to the region around E c i as shown at W = 6, and we 
find a monotonic decrease of p ex t with larger lattice sizes 
which corresponds to a localization in the whole region. 
Therefore, the destruction of the IQHE at strong mag- 
netic fields in the TBM can be well understood physically 
as due to a moving down of negative Chern integers from 
band center, caused by strong disorders, which mix with 
the positive Chern integers located at lower-energy ex- 
tended levels and eventually annihilate them, leading to 
the insulator transition of the whole system. 

Two questions remain to be answered here: Firstly, 
since only a finite-size analysis has been done so far, one 
may wonder if the above picture of the insulator transi- 
tion is still to be true in the thermodynamic limit, and 
secondly, whether such a picture also persists continu- 
ously into weak magnetic field case which is more rele- 
vant to the experimental situation. To answer the first 
question, we need to determine the thermodynamic local- 
ization length by a standard finite-size scaling method. 
[p|,p"c| The results for the same field strength as in Figs. 
1 and 2 are plotted in Fig. 3 for the case of W = 4 and 
5. At W — 4, the thermodynamic localization length 
follows a scaling behavior £ ~ \E — E c i\~ x with the ex- 
ponent x = 2.4±0.1 on the two sides of the critical energy 
E c i , which verifies the existence of a delocalization fixed- 
point shown in the inset of Fig. 1 (the extended-state 
position E c \ = —3.40 also coincides with the one pre- 
viously determined by the Chern number calculations). 
The localization length £ exhibits an another divergence 
near E C 2 ~ —2.52 (p ex t in the insert of Fig. 2 also shows a 
second peak here). Above E C 2 one finds the localization 
length to be always finite up to the band center, indi- 
cating a localized region with zero Hall conductance as 
the negative Chern integers move down to E C 2. Thus, be- 
tween E c i and E C 2 there must remain a Hall conductance 
plateau as shown in Fig. 1. With W increasing from 4 
to 5, E c i is only slightly reduced (from —3.40 to —3.48) 
whereas E C 2 moves down much more quickly (from —2.52 
to —3.0), such that the v = 1 IQHE plateau between E c i 
and E C 2 becomes narrower, and eventually these last two 
extended levels merge and disappear at a critical W c . 
Correspondingly the v — 1 plateau vanishes in Fig. 1 be- 
yond this point. Thus, in contrast to the finite-size local- 
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ization length calculations in Ref. |4[ where all the peaks 
of the localization length remain unchanged before their 
disappearance, the present thermodynamic localization- 
length data clearly indicate a moving-down and mixing 
of extended levels before they disappear, consistent with 
the Chern number interpretation for the transition. 

The magnetic flux strength <j) = 2-k/M at M = 8 - 24 
used in the previous calculations are relatively high fields. 
In order to answer the second question about the weak 
field limit, we have to reduce the field strength substan- 
tially. In the following we inspect this case by examining 
the critical disorder W c . Note that by increasing W, 
we always find that higher-energy extended levels dis- 
appear first, while the lowest one vanishes in the last 
without floating up in energy (which actually slightly 
shifts downwards). Then by following the evolution of 
the lowest extended level which can be determined from 
the thermodynamic localization length, one can decide 
a critical W c , beyond which no extended states exist in 
the whole energy regime. Wc's obtained in this way are 
shown in Fig. 4 versus 1/M, the magnetic flux strength 
in unit of 2tt. We have been able to reach a weak field 
as low as = 27r/384 or M = 384, and Fig. 4 shows 
that W c monotonically decreases with the magnetic field 
strength which is extrapolated to zero at zero magnetic 
field. This result is at variance with the conjecture made 
by Yang and Bhatt || about roughly a constant W c (~ 6) 
at M — > oo as extrapolated from the data within M < 13. 
On the contrary, W c in Fig. 4 is well fit by a law 
at M > 16 which in turn means W c oc ^/uTc (w c is the 
magnetic cyclotron energy and all quantities here are di- 
mensionless as scaled by t which is taken as the unit). If 
one estimates the Landau level broadening at weak field 
as T ~ then r oc uj c is found at the critical W c . It 
implies that once the neighboring Landau levels are cou- 
pled together due to the disorder scattering, the negative 
Chern integers near the band center will be able to flow 
all the way down to the band bottom and lead to the in- 
sulator transition. In other words, it is not necessary to 
have a strong W c comparable to the bandwidth in order 
to mix all the positive and negative Chern numbers. This 
is significant because the weak-field limit in real materi- 
als corresponds to the case that the Landau level spacing 
is much less than the bandwidth, and our results reveal 
that even in this limit the lattice effect (due to which 
one gets negative Chern numbers near the band center) 
is still crucial as far as the metal-insulator transitions are 
concerned. It suggests that the continuum model is over- 
simplified in describing such phenomena in real systems 
where the lattice structure is always present. 

Finally, we would like to make several remarks. Firstly, 
a limited floating-up picture can be still seen if the 
Landau-level filling number is fixed. The reason is that 
the number of localized states below, say, the lowest ex- 
tended level at E c \ can be enhanced by stronger disorders, 



so that when the filling number is fixed the Fermi level 
actually moves down, or, relatively the extended level 
shifts up, as already pointed out by Yang and Bhatt. || 
But such a relative "floating up" for the extended levels 
before their destruction is very limited and only within 
an order of one Landau level spacing, certainly not in- 
finite, for the lowest-energy extended level as checked 
by us up to M = 192. Secondly, as with the earlier 
work, H the present results suggest that a direct transi- 
tion from higher IQHE to a zero-Hall-effect insulator be- 
comes possible, after the IQHE plateau structure in Fig. 
1 collapses from high energy downwards to the Fermi! 
surface regime at a sufficient global phase diagram. We 
shall leave discussions about comparison with experimen- 
tal measurements elsewhere. Lastly, a further theoreti- 
cal understanding is still needed about why the negative 
Chern numbers near the band center can so easily move 
down to the band bottom once the neighboring Landau 
levels start to mix by disorder scattering. 
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Fig. 1. The Hall conductance <r# as a function of 
energy E is plotted for different disorder strength W's, 
at a lattice size 16 x 16. The inset shows the evolution 
of the v = 1 IQHE plateau at W — 4 with lattice size 
varying form 8x8 (o), 16 x 16 (•), to 24 x 24 (+). 
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Fig. 2. The density of extended states p ex t versus 
E. Disorder strength W's are the same as in Fig. 1 
with a lattice size 8x8 (about 400 random-potential 
configurations are used). The inset shows p ext at W — 4 
with lattice sizes varying as 8 x 8 (•), 16 x 16 (x), and 
24 x 24 (*). 

Fig. 3. Thermodynamic localization length £ is shown 
as a function of energy E at W — 4 and 5, respectively 
Critical behaviors near divergent points are fit by scaling 
laws (see text). 

Fig. 4. The critical strength of disorder W c versus the 
magnetic flux strength (in unit 2n) 1/M. The dashed 
curve of a \ j\[M law well fits the data at M > 16. 
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